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A non-perturbative effective model is derived for the Higgs sector of the standard model, described 
by a simple scalar theory. The renormalized couplings are determined by the derivatives of the 
Gaussian Effective Potential that are known to be the sum of infinite bubble graphs contributing 
to the vertex functions. A good agreement has been found with strong coupling lattice simulations 
when a comparison can be made. 
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The nature and eventual existence of the Higgs Boson is one of the major problems in particle physics. While direct 
searches already provide a lower bound above 100 GeV [1| for the Higgs mass, electroweak precision measuramcnts 
still suggest a light Higgs, unless new physics should be found at the TeV scale by the LHC: in that case a heavy 

r 1 Higgs would require the study of non-perturbative effects on the effective theory which describes the Higgs sector. 

; On the other hand a light Higgs is surely consistent with a perturbative treatment of the Standard Modelf SM) , but 
^p. does not rule out a large self-coupling of the Higgs sector, since a light Higgs mass has also been predicted in the 
O i' strong coupling regime by non-perturbative methods. Thus, even for a light Higgs, some important non-perturbative 
1} , effects may come into play. 

Moreover, even if the self coupling were weak, a strong Yukawa coupling is expected for the top quark, and non- 
perturbative effects cannot be ruled out in the resulting Higgs-top coupled theory. It has been recently shown[4| that 
£NJ ■ the weak couplings of the full SU(2) x J7(l)gauge theory play a very negligible role in the Higgs sector, and that the 
simple self-interacting Xtp 4 scalar theory[2|, [3| yields the same predictions of the full non-Abelian gauge theory. Thus 
we believe that the simple scalar theory, and its scalar-fermion extension^ still represent a valid starting point in the 
' study of the non-perturbative features of the Higgs sector. 

, The scalar theory has been extensively studied in the past, but few truly non-perturbative treatments have been 
• ■ reported, and most of them have failed in the attempt to give a fully consistent interacting theory independent of 
any regularization scheme. More recently the triviality of the theory has been generally accepted and nowadays the 
SM is regarded as an effective low energy model: it does not need to be weakly coupled nor renormalizable, but it 
is supposed to be only valid up to some energy scale. An energy cut-off may be used as a regulator of the diverging 
• • . integrals, and must be left finite in order to avoid to front an useless non-interacting trivial theory. In this framework 
there has been a renewed interest in the non-perturbative behaviour of scalar theories with a large but finite energy 
cut-off [f| 0|- Lattice simulations are the most reliable non-perturbative approaches: even if they fail to provide any 
5^ ■ analytical description themselves, they can be regarded as the best benchmark for testing approximate analytical 
tools. Unfortunately there are not many of such tools, and in the case of a single scalar field 1/N expansions are ruled 
out. 

The Gaussian Effective Potential (GEP)(8j-[l7| is a simple variational tool which has been often used for describ- 
ing the spontaneous symmetry breaking in the framework of the scalar theory for particle physics and condensed 
matter [l8l - [20l |. Moreover, as it has been already pointed outpll [2^|. the GEP contains information about the renor- 
malized one-particle-irreducible (1PI) n-point functions: the derivatives of the GEP are a variational estimate of the 
1PI functions at zero external momenta. Thus a direct expansion of the GEP around the broken-symmetry vacuum 
can be regarded as a an effective low-energy model with the derivatives that act as renormalized couplings. 

In this paper we study the emerging effective model and discuss the meaning of the couplings in terms of infinite 
sums of special classes of graphs. In fact the couplings are known [2lTj2o1 | to be given by the infinite sum of all the 
bubble graphs which can be drawn as perturbative corrections for the variationally optimized Gaussian Lagrangian. 
At variance with previous work the derivation of the effective model and its renormalization is carried out at a finite 
cut-off exactly, without any further approximation that would require or assume a very large cut-off which should 
eventually be sent to infinite. Thus the calculation is in the spirit of recent lattice simulations [|| 0|, and can be 
compared to such numerical calculations while providing an analytical non-perturbative low-energy effective model. 
When a direct comparison can be made, we find a very good agreement with lattice data(27j gaining confidence on 
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the reliability of the variational method. 

It is worth pointing out that, while the method is not new by itself, no previous attempt had been made to compare 
the predictions of the finite cut-off effective model with the available lattice data. For instance we derive a very simple 
analytical expression for the transition point that fits the lattice data very well up to very large couplings. We think 
that the method could be extended and used for an analytical description of the Higgs-top model which has been 
recently addressed by lattice simulation^ 

Another important issue is the existence of an intermediate mass which plays the role of a variational parameter 
and appears as an intermediate energy scale: quite smaller than the cut-off and still rather larger than the physical 
masses. Since all internal lines in bubble graphs are evaluated with such intermadiate mass, then the resulting n-point 
functions have a very weak dependence on the external momenta at low energy, and they can be regarded as constant 
renormalized couplings. In other words the physical low-energy amplitudes can be derived at tree-level by an effective 
lagrangian whose couplings already contain the sum of infinite bubble graphs. Moreover, the non-perturbative nature 
of the method does not require the self-coupling to be small. 

The paper is organized as follows: in section II the effective model is defined for a self-interacting scalar field; in 
section III the couplings are recovered as the sum of bubble expansions around the Gaussian ansatz; in section IV 
the problem of field renormalization is addressed and the nature of the critical point is studied; in section V the 
phenomenological content of the effective model is discussed; some final remarks and directions for future work are 
reported in section VI. 



II. THE VARIATIONAL EFFECTIVE MODEL 

Neglecting couplings to other fieldsQ, the Higgs sector of the Standard Model can be described by a simple scalar 
theory: in the Euclidean formalism the Lagrangian reads 

L = \d^d^+\m\<f + ^\ B <t>\ (1) 

Denoting by <p the expectation value of the scalar field ip =< (f> >, and by h the Higgs field h = <f>— <(/)>, the 
Lagrangian can be split as 

£> = &GEP ~\~ £>int (2) 

where 

£gep = \d»h d„h + ^n 2 h 2 (3) 

and Li n t = & — &gep- The GEP can be recovered (l7j as the first order effective potential for the free theory described 
by Lgep in presence of the interaction Lmt- The mass f2 of the free theory is then determined by requiring that for 
any value of the average tp the effective potential is at a minimum. A trivial calculation of the first order effective 
potential yields 

W(*0 = \™W + ^n|/ (O) + ^ 4 + ^-<p 2 Io(ty + ^ [Io(n)} 2 - l^I (il) + h(Q) (4) 
where the Euclidean integrals Iq, I\ are defined as 

Here the symbol f. means that the integrals are regularized by insertion of a cut-off A so that k < A: according to 
the well known triviality of the X(j> 4 theory the Higgs sector is regarded as an effective model with a high energy scale 
A which plays the role of a further free parameter [gj]. The GEP is given by the effective potential Vgep in Eq.(|3]) 
provided that the mass parameter Q is regarded as an implict function of ip, defined by the minimum condition (gap 
equation) BVgep I = which reads [17j 

n a = m| + ^ + ^/ (n). (7) 
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The phenomenological broken-symmetry minimum of the GEP occurs at tp = tpo where the partial derivative of the 
GEP vanishes. By insertion of the gap equation Eq.© the derivative of Vqep reads 

dVGEP 1 / 2 1 ,\ 

~w~ - 2 r " * BV ) (8) 

and it vanishes at 



fa = 



m 2 



(9) 

ip = ip 

The known phenomenology of the Standard Model requires that tpo = v = 247 GeV, and then Eq.© gives the mass 
parameter Qq = Q\ v=tpo at the minimum as a function of the bare self-coupling A^. At the minimum point tp = tpo 
the gap equation Eq.fJT]) can be satisfied by a proper choice of the free parameter m 2 B: and the theory has one only 
free parameter, i.e. the bare self-coupling (besides the cut-off A). We must stress that Vl is a simple variational 
parameter and there is no reason to believe that its value has any physical relevance. In fact the mass f2 is the mass 
of a free particle described by the unperturbed Lagrangian Cqep, while the true Higgs mass comes out from the 
1PI 2-point function T 2 which is the infinite sum of 2-point graphs that arise from the interaction Li nt . Of course a 
comparison of Eq.([9]) with the tree level perturbative Higgs mass M 2 = Xpv 2 /?> tells us that at the broken symmetry 
vacuum tp = tpo = v, and for a small bare coupling A# <C 1, the true mass is Mh ~ fio a s we expected since the 
residual interaction in Li nt becomes very small. 

In general, the exact quantum effective action T[tp] can be written as an expansion in powers of (tp — (po) around 
the vacuum expectation value ipo 



r M = i n / (v(Pi) - MP*)) 



T n (pi,.,.,p n ) (10) 



where the functional derivatives r„ are the exact 1PI n-point functions. For the theory described by the Lagrangian 
L = Cqep + <&mtj the n-point functions are the sum of infinite n-point graphs where the vertices are read from 
the interaction £j nt and the free propagator Gq(j>) = (jp 2 + CIq) arises from the optimized free-particle Lagrangian 
Lgep with Q set at the vacuum value fio according to Eq.((9]). Any physical amplitude may be evaluated as a sum 
of connected tree graphs with vertices provided by the n-point functions r„ which play the role of renormalized 
couplings (2^| . 

Now suppose, as it turns out to be the case, that the variational parameter CIq is quite large compared to the 
physical Higgs mass Mh- Then the n-point functions have a very weak dependence on the external momenta as far 
as these are small or at least pi w M^. An approximate low energy effective model can be recovered by taking these 
couplings at their zero momentum value T n r n (0) = r„(0, . . . ,0) and going back to the direct space where an 
effective Lagrangian can be written as 

L eff = ±d»h d^h + Y' 2 {Q)} h 2 - |r 3 (o)* 3 - ^r 4 (0)/i 4 + • • • (11) 

having denoted by T' 2 the sum of first and higher order contributions to the 1PI 2-point function, i.e. all the terms 
except the zeroth-order one which according to Eq.([3J is equal to — 0§ . The effective Lagrangian in Eq. tfTTj) describes 
a scalar Higgs field with a renormalized mass AfJ. = — ^(0) = il 2 — T' 2 (0) and renormalized couplings — — 1^(0), 
Xp, = — r4(0), and so on. The advantage of the effective Lagrangian in Eq. (|lip is that any low energy amplitude may 
be evaluated at tree level, neglecting all loops that have been already summed up to all orders in the renormalized 
couplings. In other words the mass Mr can be regarded as the true mass Mh of the Higgs boson, and the couplings 
gn, Xr are related to the phenomenological scattering amplitudes of the particle. 

Unfortunately we do not have the exact quantum effective action T[ip], but we can extract a variational estimate 
of the couplings r ra (0) from the GEP. In fact for constant background fields ip, ipo the effective action becomes the 
opposite of the effective potential, and if we set ip(p) — (2tt) 4 S 4: (p)p and po(p) — (27r) 4 <5 4 (p)<po in Eg. ([111)) then the 
exact effective potential reads 

^//(v) = -E^r„(o)( v -^r (12) 

n 

which is an expansion of the exact effective potential around tp = ipo. Thus the renormalized couplings can be 
calculated by the simple derivatives of the effective potential as 

r«(o) = - ■ (is) 



dtp' 



tp=<p 
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On the other hand the GEP in Eq.Q is a variational approximation to the exact effective potential V e ff and we can 
evaluate a set of approximate couplings as 



r„(o) 



d n V G 



EP 



dip'' 



(14) 



tp=<p 



Insertion in Eq. (|llj) provides a simple way to perform non-perturbative low-energy calculations by the simple eval- 
uation of connected tree graphs. In spite of the approximations, in the strong-coupling regime of the Higgs sector 
the predictions of this effective model are expected to be more reliable than perturbative calculations. In fact the 
approximate renormalized couplings in Eq. (|14p are known to be the sum to all orders of bubble graphs for the ver- 
tex functions [1^, H3) Moreover the derivatives of the GEP in Eq. ([T4"l) were shown to be a genuine variational 
approximation for the n-point functions (2lL l23l. l24| and were evaluated by several authors [U US 0] in the past. Un- 
fortunately such derivations strongly depend on the special regularization scheme and on a series of approximations 
which only make sense for a very large energy cut-off that should be eventually sent to infinite. In more recent years, 
the general consensus on the triviality of the scalar theory and the failure of any attempt to build a meaningful model 
with an infinite cut-off, has changed our view of the standard model, and the more modest aim of an effective model 
has been generally accepted as a reasonable compromise @. In this framework it would be desirable to discuss the 
variational n-point functions as defined in Eq.(|14p exactly, without any further approximation, at a given cut-off A 
which is supposed to be large but not too large, and plays the role of a physical parameter that points to the energy 
scale where new physics should become relevant. On the other hand, the resulting effective model could be easily 
compared to lattice calculations where a natural cut-off is supplied by the lattice spacing. In the next section we 
derive the explicit expressions for the variational n-point functions from the GEP, through Eq.(|14[) and by direct sum 
of the equivalent bubble expansions. 



III. BUBBLE EXPANSIONS 

The best way to calculate the approximate couplings r„(0) is by derivatives of the GEP(22[. However it is instructive 
to see that the variational approximation Eq. (ll4p is equivalent to the sum of all the tree bubble graphs that can be 
drawn for the vertex functions. By tree bubble graph we mean any 1PI graph which only contains chains of bubble 
insertions that do not make any loop: in other words the external lines of a graph become disconnected whenever 
two lines belonging to the same loop are cut. This set of graphs is a sub-class of the two particle point reducible 
graphs [29] . For instance in Fig.l the 7-loop graph (a) is a tree bubble graph while the 8-loop graph (b) is not since 
its bubble chain makes a loop. Both of them are two particle point reducible. 



(a) 




(b) 




Fig. 1: Some two particle point reducible graphs. The 7-loop graph (a) is a tree bubble graph as it becomes disconnected 
whenever two lines belonging to the same loop are cut (i.e. whenever a bubble is cut); the 8-loop graph (b) is not a tree bubble 
graph as it contains a loop of loops, and it is not disconnected by the cut of a bubble in the loop. 



In order to show the equivalence, let us write the explicit interaction L int that according to Eq.J2]) and neglecting 
constant terms becomes 

L mt {h) = [m% + + \ (m| - fig + ^) + ±<p X B h 3 + ^-h\ (15) 
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The tree bubble graphs that contribute to the 1PI 2-point function T' 2 are shown in Fig. 2 where a straight line 
represents the free propagator Go = (p 2 + ^o) 1 as derived from the definition Eq.([3]) of &gep at the vacuum ip = p . 
The first order contribution arises from the first two graphs of Fig. 2 



r 



(i) 



n 2 

"0 



A 



B 2 



A, 



■io(Po) 



(16) 



and by insertion of the gap equation Eq.© we get = which is a well known property of the GEP. Thus T' 2 is 
given by the sum of all the higher order tree bubble graphs shown in Fig. 2. 



G = ( P 2 + nl 




Xb 
' 4! ' 



r' 



-2— + -o- + -OO + <xx> 



Fig. 2: Tree bubble graphs contributing to the 1PI 2-point function. Explicit expressions for the free propagator and for the 
vertices are reported on the top. 

Let us denote by L n the n-order 1-loop graph without vertex factors, without external lines and with external 
momenta set to zero, as shown in Fig. 3. 



L 



d%p 



1 



1 



d n I 



" ! ! / a (2^) 4 ( P 2 + n\ dippy n=Uo 

and let us denote by B the bubble chain geometric expansion reported as an hatched bubble in Fig. 4 

1 



n=0 



Xb 

T L 2 



l + 4fL 2 



(17) 



(18) 



o 




= u 



Fig. 3: Graphs for the integrals L2, L3 and L4. They are n-order 1-loop graphs without vertices, without external lines and 
with external momenta set to zero. 

With the above notation the approximate 2-point coupling T' 2 is shown in Fig. 4 as a cross-hatched bubble with 
vertex and symmetry factors added at the external points, while the simple cross-hatched bubble represents the chain 
bubble sum without external vertex factors. Accordingly we can write 



/ 1 \ 3f2 2 
r 2 (0) = (3 x 3 x 2) i--\ BVo \ L 2 B = ^\ B L 2 B 

where we have made use of Eq.©. The renormalized mass Mjj = — r 2 then follows 

1 — XgL 2 



Mr = Oq — r 2 (0) = fig 



On the other hand, the derivative of Eq.([5]) yields 



d 2 V G EP fo2 1. 2 ,^ 9 2 fdn 2 X B 



and, at ip = tpo, by insertion of Eq.([2]) we get 

~ d 2 V GEP 
dp 2 



6Q 2 



dfl 2 

dp 2 



Xb 
3 



(19) 



(20) 



(21) 



(22) 



B 



1 + 



-o 



+ — GO — + 



GOO 



Lo-B 



(3x3x2)(-iA B o )%-5. 



Fig. 4: The bubble chain geometric expansion denoted by B in the main text (hatched bubble). The cross-hatched bubble 
represents the sum of all the bubble chain graphs without external vertices. The coupling F^O) is obtained by adding the 
external vertices and the correct symmetry factors to a cross-hatched bubble. 



The derivative of il 2 as a function of p 2 can be easily obtained by the gap equation Eq.© whose derivative reads 



dn 2 

dip 2 



Xb 
2 



A s / dl \ dVL 2 
~2 \citf 2 ) 'dp 2 ' 



Then the derivative of Q 2 is 



which at the vacuum ip — po becomes 



dn 2 

d^ 

dfl 2 
'dp 2 



[1 



Ab ( dl a \\ 
2 ld?F/J 



ip=ip 



2 [1 + AfLa] 



(23) 



(24) 



(25) 



Insertion of the derivative in Eq. ([22j) shows that, according to Eq. (|T4j) which defines the couplings as derivatives of the 
GEP, Eq. ([22|) yields a renormalized mass Mr which is exactly the same as that obtained by the sum of tree bubble 
graphs in Eq. ([2U|) . 

The equivalence can be extended to higher orders: the derivatives of the GEP are easily evaluated by Eq. (|2"Tj) and 
give 



9R 



d 3 V, 



GEP 



dp 3 



M 2 

2 




4^ 



'2n2 1 



d 2 n 



d^Va 



EP 



dp 4 



where the derivatives of fi 2 follow from Eq. 
Eq.dH]) 



A M « 



4\p 2 



2o2 



d 2 n 



d{p 2 ) 2 



d(<p 2 ) 2 
+ 8p 4 



d 3 n 2 



d{p 



2Vi 



(26) 



(27) 



d 2 n 2 



d(tp: 



and can be written in compact form by insertion of Eq. (|17p and 

(28) 



ip = ip 
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^B 3 L, 



d 3 n 2 



d{p 2 



ip = ip 



— Z^B-B 5 Ll — -\%B 4 L i . 



(29) 



With the above notation the renormalized couplings read 

~M 2 



Xr — X B 



9r = <^oAb 



-^t + { Vq \b) 2 B 3 L 3 



M 2 



(30) 
(31) 



(3!) +3 




(3 x 4!) 



- V (3!x3!x3!) 



Fig. 5: Tree bubble graphs contributing to the 1PI 3-point vertex function with their symmetry factors. 



r 4 = — Xr 




(4x3x4!) + 



(3! x 3! x 3! x 3!) 



(4! x 3! x 3! x 3! x 3!) 



Fig. 6: Tree bubble graphs contributing to the 1PI 4-point vertex function with their symmetry factors. 

It is not difficult to show that these couplings can be recovered by the sum of all the tree bubble graphs that can be 
drawn for the vertex functions. The sum of the corresponding geometric expansions are reported in Fig. 5 and Fig. 6. 
The tree bubble graphs that contribute to the 1PI 3-point vertex function are reported in Fig. 5 with their symmetry 
factors. The sum of the first two terms is (poXs — 3F' 2 /(po — (PqXbM r /Q,q, while the third graph yields X b (PqB 3 L 3 , 
in agreement with Eq. (f3T))) . With the same notation the tree bubble graphs contributing to the 1PI 4-point vertex 
function are reported in Fig. 6: as for the 3-point function the sum of the first two terms is XbM r /CIq while the other 
graphs can be easily shown to reproduce the right hand side of Eq. (|31[) term by term. 

The equivalence explains the content of the variational effective model and makes it even more evident that the 
GEP provides a set of genuine non-perturbative approximate couplings that represent the sum to all orders of tree 
bubble graphs for the vertex functions. 



IV. RENORMALIZED COUPLINGS AND CRITICAL POINT 



Before we can compare the predictions of the model with the phenomenology, we must address the problem of field 
renormalization. As we already said, the variational mass f2 turns out to be quite larger than the physical Higgs mass 
M/j, thus reducing any momentum dependence of the n-point functions. As a consequence we expect a negligible field 
renormalization and a physical mass M 2 which is basically equal to M% = — ^(0). However, we would like to avoid 
any further approximation, and we prefer to evaluate both the physical mass and the field renormalization factor 
Zh from the 2-point function T 2 (p) according to the usual definitions: 

r 2 (p)| p 2 = _M2 =0 (32) 



7 -i _ dr 2 ( P ) 

Thus our finding that Zh ~ 1 and Mh ~ Mr can be regarded as a check that the momentum dependence can be 
neglected, and a confirmation that the 2-point function can be taken as — [^(p)] -1 = (p 2 + Af 2 ) as it was already 
assumed without proof in previous work on the full gauge theory The same conclusion had been reached by 
numerical simulations on the lattice[6j, [7J. Accordingly, the physical renormalized couplings must be defined as 

X h = -Z 2 r 4 (0) - Z\\ R . (34) 



g h = -Z 3 h /2 T 3 (0) = Z%*g R 



3/2 



(35) 



The full 2-point T 2 (p) function may be recovered by the sum of the bubble expansion in Fig. 4 with the loop L 2 
replaced by the function L-^ijp), defined by the same graph in Fig. 3 with the external momentum restored 



LM =m 



1 



1 



I a (2^) 4 [( P + k)i + n 2 ] [fc2 + r> 2 ] 

Then according to Eq. (|19|) the 2-point function reads 

3tt 2 \ B L 2 {p) 



r 2 (p) = (^ + ng) 



2 1 



l L 2 {p) 



(36) 



(37) 



The same result can be found by the covariant Gaussian approximation[24|, or by direct functional differentiation of the 
Gaussian effective action[2lj, and can be shown to be a genuine variational bound of the exact 2-point function [2ll [23| . 
According to the definitions in Eqs. (f3"2"l) . ([3"3")) the physical Higgs mass Mh is obtained as a solution of the equation 



Ml = n% 



1 - X B L 2 (p) 
l + ^L 2 (p) 



(38) 



J p2 = _ M 2 



while the field renormalization constant is given by 



Z7 1 = 1 - ^ 2 A 



A B 



dp 2 



(P) 



(39) 



The integral in Eq. (|36p must be evaluated inside the four-dimensional hyper-sphere k < A, so that the usual 
Feynman formula cannot be used. However, by a tedious calculation the integral can be shown to yield the following 
exact result 



1 



32ttV 



O 2 4- A 2 
A 2 + P 2 log^±A 

"0 



Hp) 



where the integral I(j>) 



Hp) 



Q o +A2 dx 



— \lix-p 2 ) 2 + Vfi 2 , 

X v 



(40) 



(41) 
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is given by 

i(p)=p 2 ti(ta-ti)-log 

with the variables U,t± defined according to 



-1 \- x l+U\ log 



1 

2^ 



{t 2 -t-){t 1 -t+) 



(*2 -*+)(*!-*-) 



(42) 



(43) 



ti = 



(44) 



2pft 



a 2 - p 2 + ft 2 , + J vng + (ng + a 2 - p 2 



(45) 



In spite of the appearence, the loop function L2(p) in Eq. (140p has a logarithmic divergence only: in fact for 
Mh, A we can write, up to order 1/A 2 



L 2 (p) 




1 + log 



20 2 



(46) 



yielding at the point p 2 

^2(p)| p 2 = _ M . 



M, 2 



16tt 2 



1 + 2 log 



A 



arctan 



M hv /4r> 2 -M 2 
2Q 2 , - Ml 



(47) 



We notice that as far as A is finite, large but not too large compared to Qq, the constant terms must be retained as 
they can be of the same order as the logarithm log(A/Oo)- From Eq.(|47|l we see that L2(iMh) is analytical in the 
limit Mh — > and its limit value 1/2(0) is 



£ 2 (0) = 



1 



16tt 2 









- 1 







(48) 



which is exactly the same large A behaviour which would come out from a direct calculation through the definitions 
of L n and J by Eqs. (fT7j).(|5|>. 

The limit Mh — > is very important as it turns out to be the critical point of the Higgs sector, since the vanishing 
of the mass may be associated to a change of sign for the second derivative of the effective potential at its stationary 
point which becomes a maximum. The exact equation for Mh, Eq . (|3"8"|) , can be easily seen to admit the solution 
Mh = whenever the condition 



Ar = 



1 



L 2 (0) 



(49) 



is fulfilled. Thus we observe that, for a large but finite cut-off, the bare self-coupling Xg does not need to be small at 
the critical point: assuming that fi / A is small enough we can make use of Eq. (|48[) and write the critical condition 
Eq.® as 



Q = -7= e" 



(50) 



We can check that for a coupling as large as Xb ~ 871" 2 « 80 we still find flo/A w 0.22 which is small enough to be 
consistent with the use of the approximate Eq . (|48jl for £2(0). 

From Eq. (|50p we see that there is no way to keep fio finite when the cut-off A is sent to infinity unless the 
bare coupling As is taken to be infinitesimal. This is just what happens in the autonomous renormalization(l7l| 
where the bare self-coupling is taken to be = 1/^2(0) while £2(0) diverges logarithmically. That is exactly the 
required condition Eq. (l49p for the vanishing of the Higgs mass. Thus it is not a great surprise that in autonomous 
renormalization the GEP predicts a vanishing Higgs mass[21]. 
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From Eq. ([50|) we also see that the GEP predicts a weak first order transition because Oo does not vanish at the 
critical point. According to Eq.© at the minimum of the effective potential the vacuum expectation value of the 
scalar field is <pQ = 30 2 ,/ Xb, and if Oo does not vanish at the transition, the expectation value ipo jumps to zero at 
the critical point. That seems to be an unavoidable shortcoming of the GEP, since the transition is believed to be 
continuous. However by Eq. (|50[) we see that the jump Aip of the expectation value at the transition is small 

v^A 

Aip = - e x b . (51) 
\fe\ B 

It reaches its maximum at Xb = 1&tt 2 ~ 158 where Aip/ A pa 0.05 and it goes to zero in both the limits Xb — > and 
Xb — > oo. Thus we may neglect the jump in most cases, as far as we do not go too close to the transition point. That 
is not a major problem at all, as we know that the physical range of interest cannot be too close to the critical point 
where the Higgs mass eventually vanishes. We cannot use the GEP at criticality, but we expect that the order of the 
transition should not affect the behaviour of the resulting effective model provided that we do not reach the critical 
point. Actually we will see in the next section that the vacuum expectation value ipo seems to be continuous up to 
the transition point for any reasonable plot resolution. 

Whenever A ^ Oo we can use the approximate Eq. (J^HJ) f° r L 2 (p) hi order to get the field renormalization constant 
Z h . Insertion of Eq.(|3S|) in Eq.® yields 



1 



X B (2 + x 2 ) 2 



96tt 2 (4 
where x = Mh/^o an d 



1 - 



32tt 2 f(x) 



(52) 



f(x) = L 2 (iM h ) - L 2 (0) 



16tt 2 



2 - 



V4 - x 2 



■ arctan 



(53) 



We can simply check that Zh pa 1 in the small-coupling limit and close to the critical point. That is quite obvious in 
the small-coupling limit Xb — > where Mh ~ Oo and x — > 1. In this limit Z^ 1 — 1 pa 2 • 10~ 3 • Xb, which is vanishing 
small. In the opposite strong coupling limit we can explore the critical range where Mh — > (x — > 0) and, according 



to Eq.gU), X B L 2 pa 1. In this limit 96n 2 f(x) 



and Z,; 1 - 1 sa A b /(144tt 2 ) pa 0.7 ■ 10" 3 • Xb- This is just one per 



cent for a coupling as large as X B sa 15. In actual calculations we have never found values larger than few per cent in 
the broken symmetry phase, in agreement with previous numerical findings |(| @|- 

Quite interesting, around criticality we find a hierarchy of energy scales with an Higgs mass Mh quite smaller than 
the mass parameter Oo which is supposed to be much smaller than the energy cut-off A. The physical meaning of the 
intermediate scale Oo is related to the vacuum expectation value of the scalar field v = 247 GeV according to Eq.© 
which reads Oq 



Vy/Xs/3- For a strong coupling Xb 
we like by a large cut-off. 



20 we get Oq pa 640 GeV while Mh can be made as small as 



V. RENORMALIZED COUPLINGS AND PHENOMENOLOGY 



In order to make contact with Monte Carlo calculations [6J, [jj we take all energies in units of the cut-off A. In these 
units the vacuum expectation value ipo = v/A can be seen as representing the inverse of the cut-off in units of v = 247 
GeV. In fact at the critical point the vanishing of ip Q can be regarded as a restoration of symmetry at a fixed cut-off, 
or as the effect of an infinite cut-off on a fixed vacuum expectation value. Following Ref.|6] we take a constant value 
for the bare self-coupling Xb and change the bare mass m 2 B up to the critical point which is reached when Mh vanishes 
at mr B = m 2 . In the broken-symmetry phase we measure the distance from the critical point by the adimensional 
parameter r = |1 — m 2 B /rn 2 \. The physical Higgs mass Mh and the vacuum expectation value ipo are reported in Fig. 7 
up to the transition point for a moderate Xb — 10. We observe that the first-order jump of ipo is so small at the 
transition that it cannot be even seen in the plot at any reasonable scale. 

A direct comparison with numerical simulations requires that the correct scale factor should be fixed between 
energies. In lattice calculations energies are usually taken in units of the inverse of the lattice spacing a which also 
provides a natural cut-off. Thus we expect that the scale factor c = Aa should be of order unity, but we have no 
direct and unique way to determine it. For instance we could require that the unit cell of the inverse lattice has the 
same extension of our four-dimensional hyper-sphere p < A, and obtain c = 2 5 / 4 y / 7r pa 4.2. Here we prefer to make a 
more empirical choice and take the scale factor c as the ratio between the critical bare masses m C} assuming that the 
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Fig. 7: The physical Higgs mass Mh (solid line) and the vacuum expectation value ipo (dashed line) as functions of the 
adimensional parameter r = |1 — m B /m^.\, in the broken symmetry phase, for a moderately strong bare coupling Xb = 10. 
Data points are the lattice simulations of Ref.[27l| for the symmetric phase, reported assuming Aa — 4.38. 



critical masses should be the same if the energy scales are correctly handled. For a moderate Xb ~ 10 — 100 there are 
not too many data to compare with: some numerical results have been reported for the symmetric phase in Ref. (27| , 
and a comparison of the critical bare masses for Xb = 10 yields the empiric scale factor c = 4.38. 

Once the scale factor has been fixed we may compare our results with the lattice data. The critical point m 2 c is 
shown as a function of A^ in Fig. 8 where the lattice data of Ref.[27j have also been included for comparison. Our 
result fits the lattice data very well up to Xb ~ 150. For larger couplings the GEP underestimates the strength of 
the critical m B as a natural consequence of the predicted first order transition. On the other hand we may extract 
from the GEP a yet simpler analytical approximate function that interpolates the data very well for a moderate bare 
coupling. In fact for a moderate coupling Xb < 100 we already know that, according to Eq. (|50j) . the ratio Clo/A is 
small enough at the critical point. Therefore we can neglect powers higher than il§/A 2 in the expansions of £2(0) 
and I (fl ), and inserting Eqs. (|5D|) . © in the gap equation Eq.([7|) we gain the simple result 



2 

m n — 



32ir 2 



1 — e 



(54) 



As shown in Fig. 8, the approximate Eq. (|54p gives the correct critical point up to A^ ~ 200, while for stronger couplings 
the strength of the critical m B is overestimated. We observe that Eq. ([5^|) has an essential singularity at the point 
Xb = and thus it cannot be recovered by any perturbative expansion in powers of Xb- Despite being very simple, 
the analytical result in Eq. ([M|) is a genuine non-perturbative prediction, and provides an example of the capabilities 
of the variational method. 

With the same scale factor, the available lattice data[27j for the mass have been inserted in Fig. 7 for comparison, 
and again they are in good agreement with the present variational calculation. Both calculations are consistent with 
a square root behaviour Mh ~ y/r. Furthermore, once the correct scale factor has been inserted, the numerical values 
seem to be almost the same, with our data slightly larger than the lattice predictions. That was not entirely expected 
as in the lattice simulation of Ref.[27| the renormalized mass was measured in the symmetric phase while our results 
are for the broken symmetry phase, and the critical behaviour is expected to agree in the two phases up to constant 
factors: integration constants that are usually different in the two phases. Some perturbative arguments have been 
reported in the past giving evidence for the equivalence of the integration constants in the two phases[30], while lattice 
data have shownQ that, for a very strong coupling Xb = 600, the constant factors are different, with the broken 
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Fig. 8: The critical parameter m 2 B = as a function of the bare coupling As. The solid line is the result of the present 
variational method. The squares are the numerical lattice data of Ref.[27| scaled according to Aa — 4.38. The dahed line is the 
outcome of the approximate analytical function reported in Eg . (154[) . 



symmetry mass being almost twice the mass of the symmetric phase. Thus we argue that the slight difference in Fig. 7 
can be taken as a measure of non-perturbative effects, still small for Xb — 10. Actually we have checked that when 
Xb increases the difference also increases with the effect becoming quite large for Xb = 50 already. Unfortunately, 
as discussed at the end of the previous section, we cannot allow Xb to become too large at the transition point since 
the transition is predicted to be first order by the GEP. That precludes a direct comparison with the lattice data of 
Ref.[7] in the broken symmetry phase. 

In Fig. 9 the physical adimensional renormalized couplings 3M 2 /v 2 , gh/v and A^ are reported. For a large r (i.e. a 
small cut-off) they all tend to the same bare value Xb- At the critical point r — > they all seem to vanish yielding a 
trivial model. However it should be observed that the renormalized coupling Xh is not equal to 3M%/v 2 as it would 
be suggested by the standard tree-level relation. Such a condition can only be satisfied if ss f2 in Eq. ([§]), i.e. 
for a very small bare self-coupling or an unphysically small cut-off (large r). The same adimensional renormalized 
couplings are reported in Fig. 10 for the critical range r < 1. For comparison in the same Fig. 10 some data points 
from Ref.[27| are reported as an indirect measure of 3M 2 /v 2 on the lattice. According to scaling, these points can be 
assumed to be equal to the value of A_r in the symmetric phase. In fact for A_r we could not find any direct simulation 
data in the broken symmetry phase, nor any direct lattice measure of the renormalized couplings Xh and gh/v. 

For small values of the parameter r (i.e. for a large cut-off) all the couplings decrease and show a logarithmic 
behaviour. We observe that A^ is always larger than its perturbative value iM 2 /v 2 , yielding a strongly coupled 
effective model even when the Higgs mass falls well below the electrowek energy scale v. Moreover the 3-point 
adimensional coupling gh/v and the 4-point vertex Xh turn out to be different at variance with the case of a weak 
coupling where they are the same. For a very light Higgs (close to the critical point), we may assume the autonomous 
condition Eq. (|49p . take B w 2/3 in Eqs. (|3"0"|) , (|3ip and for a large cut-off write the ratio gh/{Xhv) as 



vX h 



' (55) 



8tt 2 



We observe that the ratio between the couplings is gh/(vXh) ~ 0.25 for any moderate bare coupling, and still smaller 
if As gets very large, to be compared with the weak coupling limit gh/(vXh) ~ 1- The weakening of the 3-point 
coupling would give less chances of finding bound states in the Hig gs s ector, since it is the 3-point vertex that can 
give a binding contribution in a strongly self-coupled Higgs sector [2L l3l|. 
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Fig. 9: The physical adimensional renormalized couplings 3M%/v 2 (lower solid curve), gu/v (central dashed curve) and \u 
(upper dotted curve) as functions of the adimensional parameter r = |1 — m%/m^\, in the broken symmetry phase, for a 
moderately strong bare coupling Ag = 10. The couplings merge in the large r limit (small cut-off). 



VI. DISCUSSION 



The effective model that emerges variationally from the GEP could provide a reliable way to describe the experi- 
mental data, together with numerical lattice simulations. In fact we have shown that the non-perturbative predictions 
of the GEP can be as effective as numerical simulations, provided that the comparison is not pushed too close to the 
critical point where the GEP is known to fail. That precludes the study of extremely strong couplings Xb ~ 100—1000 
as according to Eq. (|5ip the jump of the vacuum expectation value tp would not be negligible at the transition point. 
Unfortunately that has not allowed us to compare our predictions with some numerical data@ on the broken-symmetry 
phase that have been reported for As = 600. 

On the other hand most of the numerical simulations that have been reported for moderate couplings An ~ 10— 100 
happen to be too far from the relevant physical range of parameters as the cut-off is usually too small [23, |32|- That 
represents a common problem of numerical simulations since the finite size of the sample limits the accuracy when 
the correlation length £ reaches the size L of the sample. Usually L = na where a is the lattice constant and n is a 
small integer of order ten at most, therefore taking Mh = l/£ and A ss 1/owe get a bound A < nMh- For a light 
Higgs Mh ~ 100 GeV and for n — 10 we get a cut-off A < 1 TeV, too small for exploring the Higgs sector. Actually, 
by a variational argument, a threshold has been predicted at A « 3.5 TeV for the existence of a strongly interacting 
light Higgs[|. 

However some simulation data have been reported for the symmetric phase and a moderate strong coupling[27], and 
have been found in perfect agreement with the predictions of the present method. Larger lattices should be studied 
for a comparison in the physical interesting range of parameters, while the real experiments seem to be an interesting 
and viable way to test the predictions of the present effective model. 

It must be mentioned that before a full comparison can be made between the real phenomenolgy and the effective 
model, the couplings to fermions and gauge fields must be included in the simple scalar theory. Inclusion of gauge 
fields is not a major issue as we have already derived the GEP for the full SU(2) x U(l) non-Abelian gauge theoryjj], 
and shown that the weak gauge couplings do not play any relevant role in the Higgs sector. The couplings with 
fermions are expected to be more important and, since they are proportional to masses, the Top quark should be 
considered at least. The present variational method can be extended to the Higgs- Top model that has been recently 
studied by the GEP [33] and by numerical lattice simulation[6|. While we do not expect any dramatic change of the 
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Fig. 10: The physical adimensional renormalized couplings 3M 2 /v 2 (lower solid curve), gn/v (central dashed curve) and \h 
(upper dotted curve) as functions of the adimensional parameter r = |1 — m%/m%\, in the critical range of the broken symmetry 
phase, for a moderately strong bare coupling Ag = 10. Note the logarithmic decreasing of the couplings for r — ► 0. Data points 
are indirect lattice simulations for 3M%/v 2 , and have been obtained by scaling from the symmetric phase data of Ref . [27l| . 



general phenomenology with respect to the present simple theory, the resulting effective model should give a more 
reliable description of the Higgs sector. 
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